Diamagnetism and flux creep in bilayer exciton superfluids 
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We discuss the diamagnetism induced in an isolated quantum Hall bilayer with total filling factor 
i^T = 1 by an in-plane magnetic field. This is a signature of counterflow superfluidity in these 
systems. We calculate magnetically induced currents in the presence of pinned vortices nucleated 
by charge disorder, and predict a history-dependent diamagnetism that could persist on laboratory 
timescales. For current samples we find that the maximum in-plane moment is small, but with 
stronger tunneling the moments would be measurable using torque magnetometry. Such experiments 
would allow the persistent currents of a counterflow superfluid to be observed in an electrically 
isolated bilayer. 
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I. INTRODUCTION 

Superfluiditjff' is a spectacular form of hydrodynamics 
involving dissipationless flow, metastable circulation, and 
quantization of circulation. It occurs in liquid helium and 
cold atomic gases, where it is associated with the con- 
densation of many bosons into a single quantum state. 
Whether condensates of quasiparticles, such as excitons, 
would also be superfluids has been discussed for many 
years, with debate over the physical manifestationa^ltSl 
of superfluid hydrodynamics for quasiparticles, the (re- 
lated) role of symmetry-breaking perturbations, and the 
significance of interactions and thermal equilibrium^j. 

This question can now be addressed experimentafly, 
with emerging evidence for the condensation of quasi- 
particles including excitons, polaritons, and magnons. A 
particularly interesting system is the quantum Hall bi- 
layer at total filling factor 1/^ = 1 (see Fig. [ij. This 
consists of two closely spaced quantum wells, each con- 
taining a two-dimensional electron gas, subjected to a 
strong perpendicular field B± so that the lowest Landau 
level in each layer is half-filled. The tunneling between 
the two layers is weak compared to the Coulomb energy 
scale. When the interlayer separation d is of the order of 
the magnetic length = {h/eBj_y^^ , the ground state 
is believed to be a Bose-Einstein condensate of interlayer 
excitonJi^HIII caused by the Coulomb attraction of elec- 
trons and holes across the layers. (Note that the holes are 
unfilled electron states of the lowest Landau level of the 
conduction band.) Flows of these condensed excitons cor- 
respond to dissipationless counterflowing electrical cur- 
rents in the two layer^^^. Initial evidence of dissipation- 
less transport came from interlayer tunnelingj ^^ l ^^l which 
exhibits non -zero interlayer currents at negligible inter- 
layer volt age j- ^ ^ I ^ s^ o\ xhis regime persists up to a critical 
current, above which dissipation increases rapidly. 

The interpretation of transport measurements as 
evidenc^^U for exciton superfluidity is complicated by 
parallel charge tr ansport channels, the injection and re- 
moval of electrongZESESl contacts, as well as pos- 



sible dissipation in the leads. Here, we show that mag- 
netometry on isolated bilayers could provide direct ev- 
idence for exciton superfluidity, without the complica- 
tions inherent in transport studies. We predict that, at 
low temperatures, the bilayer shows a history-dependent 
susceptibility. Changes in the in-plane field lead to a 
persistent diamagnetic moment that is, however, not in- 
duced by fields present when the condensate forms. This 
discrepancy corresponds to that between the moments 
of inertia of normal and superfiuid helium inferred in a 
torsional oscillator experiment. 
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FIG. 1. Schematic diagram of a bilayer of electrons in the 
quantum Hall regime due to a strong perpendicular field B±. 
Applying a weak in-plane field B causes counterflowing charge 
currents as a diamagnetic response. For an isolated bilayer the 
current loop closes at the edges of the sample by interlayer 
tunneling. In an excitonic picture, this corresponds to the 
flow of neutral excitons along the bilayer with recombination 
at the edges. 

The in-plane magnetic susceptibility of the exciton 
condensate has previously been considerecP'SiHSei jjj dgan 
bilayers. However, there are very good reasons to 



expecl! ^'"'^" ' that disorder is essential for understanding 
the observed transport properties of the bilayer at low 
temperatures. In particular, for quantum Hall bilayers, 
charge disorder nucleates vortices in the condensate^SESl 
Therefore, we will study in this work the magnetic prop- 
erties of a quantum Hall bilayer with disorder-induced 
vortices. We first consider the zero-temperature limit, 
and show how dissipationless supercurrents lead to long- 
lived diamagnetic moments. In current samples, the re- 
sulting moments are small [Eq. (16)], but for samples 
with stronger tunneling the moments can be compara- 
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ble to the Landau diamagnetism. We also consider the 
effects of non-zero temperatures, and estimate the dissi- 
pation due to thermal creepP^ of the in-plane flux. This 
dissipative mechanism gives rise to a non-vanishing re- 
sistance for the counterflow supercurrents, and causes 
the diamagnetic moments to decay in time. Based on 
the parameters of current experiments, we expect dis- 
sipation to be significant, which is consistent with ex- 
perimental resistance measurements. The fiux motion 
in this regime can be studied in transport experiments. 
Samples deeper in the condensed phase, however, would 
show persistent counterflow supercurrents on laboratory 
timescales, whose presence and eventual decay could be 
observed by magnetometry. 

The remainder of this paper is structured as follows. In 
Sec. [nj we outline the model we consider for the behav- 
ior of long-wavelength counterflow supercurrents in the 
disordered bilayer. In Sec. |III[ we apply this model to 
calculate the in-plane susceptibility of the bilayer in the 
zero-temperature limit. In Sec. IV we discuss the ex- 
tension to finite but low temperatures. Finally, in Sees. 
[V| and |VI[ we present numerical estimates for the mo- 
ments and resistance in experiments, and summarize our 
conclusions. 



II. MODEL AND BACKGROUND 

In this work, we consider a quantum Hall bilayer at to- 
tal filling factor z^t = 1, in the presence of the disordered 
electrostatic potential originating from the dopants, set 
back at a distance dj,- For definiteness, we adopt the 
"coherence network" picture of Fertig and Murthy,!^ al- 
though our theory here is more general. In this picture, 
the random Coulomb field from the dopants creates pud- 
dles of normal electron liquid surrounded by channels 
of excitonic superfluid. In quantum Hall ferromagnets, 
physical and topological charge are related so that the 
charge puddles nucleate vortices, known as merons, in 
the superfluid channel^23MI, We estimateci^ that typi- 
cal disorder strengths induce on the order of one unpaired 
vortex per puddle so that the correlation length of this 
disorder is ^ « ~ 100 nm. 

The energy associated with the excitonic supercurrents 
in the channels is, in terms of a superfluid phase 77,---- 



|V77 + qp -tcos(77-|-6'o) 



(1) 



where the first term is the kinetic energy of the coun- 
terflow supercurrent with superfluid stiffness ps, and the 
second term is the energy of the interlayer tunneling cur- 
rents with tunneling strength t. The merons introduced 
by the charge disorder give rise to the random field 9q. 
For the coherence network, both ps and t should be renor- 
malized by the area fraction of the superfluid channels. 

Eq. (0 is written in a gauge where the vector potential 
of the in-plane field B is zero perpendicular to the layers 



and non-zero parallel to the layers. This field t hen ap- 
pears in the kinetic energy, inducing a wavevect o r^ * 1 2 1 3 5 ] 



q ^ {B X z)ed/h, 



(2) 



where d is the interlayer separation. The counterflow and 
tunneling currents are seen to be 
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(V77 + q)= jp+jd, 



it = —sm{r] + eo), 



(3) 
(4) 



so the magnetic held induces a diamagnetic contribution 
jd — epsQ/h to the counterflow supercurrent jd, but does 
not appear explicitly in the coherent tunneling current jt- 
We can make a gauge transformation so that the vector 
potential is non-zero only perpendicular to the bilayer, 
in which case the phase transforms to rj ^ rj + q ■ r = 9: 



t cos{9 — q ■ r + 6q) 



d^r 



(5) 



This form for the energy functional highlights the fact 
that the in-plane magnetic response of the bilayer van- 
ishes if there are no current loops that encircle the field — 
the field disappears from Eq. ^ when t — 0. (A non-zero 
diamagnetic susceptibility remains possible in multiply- 
connected geometrieJ^ or at non-zero frequencj^^. We 
comment on the extension of our results to the Corbino 
disk geometry in Secfv}) 

At zero temperature, the state of the bilayer is deter- 
mined by minimizing Ec^. Q, which gives 



X^jV^r, + sin(ry + ^o) = , Xj = {ps/t) 



1/2 



(6) 



When ^0 = this is the pendulum equation, containing 
the Josephson length Xj: the characteristic lengthscale 
over which excitonic supercurrents decay by interlayer 
tunneling. This scale is estimated to be on the order 
of a few microns in experiments. We will discuss below 
the decay of excitonic supercurrents by tunneling in the 
disordered case, which involves a different lengthscale. 

It is useful to recall, for comparison withour treat- 
ment of the disordered case, the responsJ^'^^'^^®^ of 
the clean model {9q = 0) to an in-plane field. For 
small fields, the ground state will minimize the tun- 
neling energy, and so 9 « q ■ r or rj sa 0. This is 
the commensurate state, in which the field induces a 
long- wavelength counterfiow supercurrent jcf ~ Jdi as in 
Eq. ([3|. Thus there is an in-plane magnetic moment of 
M\\ — jddLxLy ~ dLxLyXoB / iJ,o, with the susceptibility 



Xo = Pajd/B = p.oe^psd/r? 



(7) 



We see that the diamagnetic moment M\\ = xoB in- 
creases linearly with the in-plane field B in the commen- 
surate state. 

In an isolated bilayer, the counterfiow currents must 
vanish at the ends of the sample. This occurs in the 
commensurate phase because, as dictated by Eq. (|6|, 
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the diamagnetic currents in the bulk of the sample are 
eliminated by tunneling over a region of size Xj near the 
sample boundaries (Fig. [T]). However, since the phase 77 
is a periodic variable the maximum current which can 
recombine in this way is given by eps\^ri\/h ^ eps/hXj. 
If the diamagnetic current present in the bulk exceeds 
this value, then a net winding of the phase 77, in the form 
of Josephson vortices, enters from the boundaries. This 
occurs at the field 



Bq ^ h/edX,j 



(8) 



Above Bq the system is in the incommensurate phase in 
which the kinetic energy of the counterflow supercurrents 
dominates so that 0K,QoirjK,—q-r. The net winding 
in T] along the sample implies that the counterflow su- 
percurrents have weak oscillations around zero along the 
sample, and the diamagnetic moment is small. The field 
-Bq marks the field above which the diamagnetic suscepti- 
bility rapidly decreases from xo ^ Josephson vortices fill 
the system and compensate the diamagnetic contribution 
in Eq. ([3j. 

While in the clean limit the phase twists nucleated at 
the boundary propagate into the bulk, in the disordered 
case any such phase twists can be pinned by the disor- 
der. They can thus be kept out of the bulk of the sample, 
which can continue to contribute the full diamagnetic 
susceptibility xo- To describe this effect quantitatively 
we shall use t he co llective pinning theory we have pre- 
viously appliecP^'^SI transport experiments in zero 
field (g = 0). 

The starting point for understanding the transport ex- 
periments is to note that Eq. ([T]) is a random-field XY 
model. The energy HcS consists of the competition be- 
tween the tunneling energy, which is minimized by a spa- 
tially varying superfluid phase 77 over the disorder corre- 
lation length ^, and the superfluid stiffness, which is min- 
imized by a uniform superfluid phase. For the bilayer, the 
disorder correlation length is much shorter than the clean 
tunneling length: ^ ^ Aj. In this regime, the superfluid 
stiffness dominates at short scales, up to the Imry-Ma or 
pinning length 



La 



(9) 



where the two energies balance. We est imat^^ll that w 
10 - 100 /im and « 20 - 100 mK in experiments. 
Beyond this scale, the phase rotates to take advantage 
of the tunneling energy. Therefore, we can interpret the 
ground state as consisting of randomly polarized domains 
of size id- The total tunneling current in each domain 
is zero in the ground state because the random field 
gives a current density jt [Eq. Q] of random sign within 
each domain of constant 77. 

However, if we drive the system by injecting currents 
phase twists enter from the contacts, leading to configu- 
rations with non-zero tunneling currents. Injected coun- 
terflow current decays into the bulk via tunneling, i.e., 
recombination of excitons. This tunneling current is sup- 



plied by rotating the phase of the domains near the con- 
tacts, leaving domains in the bulk in their ground state. 
The maximum coherent tunneling current that can be 
supported by each domain is given bj^^Il 



= eps/h. 



(10) 



If all the domains near the contact are rotated to supply 
this maximum tunneling, then we see a uniform tunnel- 
ing current jt until all the injected counterflow current 
has decayed by tunneling. Then both the tunneling and 
counterflow currents are zero in the bulk. This can be 
described by the continuity equation for the current: 



(11) 



where L\ is the size of the domain and /j is the tunneling 
current across the layers in the domain. The sign of /d is 
determined by what is necessary to reduce the counter- 
flow current. In a one-dimensional geometry, this gives a 
counterflow current profile that decays linearly into the 
bulk from the edge. For instance, for an injected current 
of —ji at a; = 0, the current profile in the cc-direction is 



-ji + hx/Ll 







Q < X < xo 

X > Xq . 



(12) 



where the point xq = jiL^/Id marks the boundary of the 
region from the sample edge where coherent tunneling 
occurs to reduce the counterfiow current. This current 
profile is a critical state similar to the Bean critical state 
in superconductor^^. In this state, there is a region near 
the contact over which the density gradient in the soliton 
train (introduced by the injected current) balances the 
pinning force arising from the tunneling. 



III. SUSCEPTIBILITY WITH DISORDER: 
ZERO-TEMPERATURE LIMIT 

We now consider the response of the disordered sys- 
tem to an in-plane field at zero temperature. We will 
see that the disordered system has a different diamag- 
netic response from the clean system. It does not have a 
commensurate-incommensurate transition controlled by 
an intrinsic lengthscale Xj [see Eq. (8|]. Instead, we find 
a saturation phenomenon for the diamagnetic moment 
at a field Be [sec Eq. ( [T5| ] controlled by the properties 
[Eq. ^ and ( [lO| )] of the phase-polarized domains, i.e., 
the disorder-pinning of the applied flux. This is followed 
by a depinning of the polarized domains and a strong 
suppression of the diamagnetic response at a higher held 
Bc2 > Bq [Eq. ( 17)] controlled by the disorder correlation 



length ^. Such behavior is similar to that of the mixed- 
state of a superconductor, but very different from that 
previously predicted for exciton condensates. 

We require the response of the superfluid phase 77 to an 
in-plane fleld B, corresponding to a non-zero q in Eq. (0. 
To obtain this, we note that Eq. ([6| for the phase rj does 
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not depend on the in-plane field. Alternatively, we can 
see that the superfluid phase is coupled to the field in 
Eq. ([l]) as an integral over q.^rj, which can be written 
as a boundary term for rj. Thus the current distribution 
in a field can be related to that in zero field with shifted 
boundary conditions. Thi s obse rvation allows us to ap- 
ply the critical-state model-^'-^' in zero field, as reviewed 
above, to calculate the current profile in non-zero field. 

For definiteness, we consider a rectangular sample with 
dimensions and Ly, with the in-plane field in the y- 
direction. Thus, counterflow currents jci{x) flow in the x- 
direction, and tunneling currents jt{x) in the z-direction 
(see Fig. [T]). For an isolated bilayer, jcf{x) = at a: = 
and Lx. We can decompose jd into paramagnetic and 
diamagnetic parts as in Eq. ([s]). The paramagnetic part, 
jp = eps^Tj/h, obeys the boundary condition ^p.^rj^^^o = 
j-p.x\x=L^ — —jd- Since the phase rj obeys the same 
equation, Eq. ^ , in the bulk irrespective of an in-plane 
field, jp has the same profile as in a system in zero field 
(when jcf — jp) with a current of —jd flowing across 
the boundaries. The full counterflow current jd for the 
isolated bilayer in a field is recovered simply by adding a 
uniform +jd to this zero- field profile oi jp(x). 
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FIG. 2. Profiles of the counterflow (solid lines and circles, left 
axis) and tunneling (dashed, right axis) currents induced by 
applying an in-plane field _B to a rectangular bilayer (Fig.[l]), 
at zero temperature. The lines show results from the critical- 
state model, and the circles show numerical results from the 
microscopic Eq. ^ . The top (bottom) panel shows the profile 
below (above) the critical field [Eq. (14l]. The tunneling 
currents have opposite signs in the two halves of the sam- 
ple (not shown). The simulation (circles) is a 200-site lattice 
model in one dimension with ta^ / ps — 0.4, where a is the lat- 
tice spacing and ga = 6 (top) and 12 (bottom). The disorder 
is uncorrelated from site to site (^ « a). 

In the critical-state model as discussed above, the 



linearly in space from —jd at the edge to zero. In other 
words, the profile of the paramagnetic current jp{x) near 
the X — Q edge is given by Eq. ( 12 1 with ji = jd- The dis- 



tance xa which gives the width of the region where coher- 
ent tunneling occurs is determined by the distance where 
the total counterflow current vanishes: jpixo) + jd = 0. 
This criterion gives a distance of 



xo = jdLl/Id = qLl 



(13) 



from the edge of the sample. The total current profile is 
shown in Fig. [2] (top). The fi gure also shows the tunnel- 
ing current density \jt\ — Id/L\ which is at its maximal 
value up to the edge of the tunneling region at qL\ from 
(either) edge of the sample. Beyond zq, there is no 
paramagnetic contribution to the counterflow current. 

This current profile is similar to the one depicted in 
Fig. [T] A diamagnetic counterflow is generated by the 
in-plane field in the bulk of the bilayer. However, an 
isolated bilayer must have zero counterflow current at 
the edges. In our critical-state model, this is achieved 
by coherent tunneling via phase-polarized domains. The 
width of this tunneling region is determined by the size of 
the domain and the coherent tunneling current that 
each could support I^- We see that a stronger fleld gives 
a higher diamagnetic current and so more domains must 
be involved. In other words, the size of the tunneling 
region increases with the in-plane field, as can be seen 
from Eq. (13 1: xq on q cc B. This should be contrasted 



with the decay length Aj for the commensurate state of 
the clean model which is an intrinsic scale that does not 
vary with the field. 

As the parallel field B is increased, the diamagnetic 
current in the bulk increases linearly with B. This re- 
quires a wider tunneling region so that all the current 
can decay to zero at the edge. So, the tunneling re- 
gion increasingly penetrates the bulk. This penetration 
is complete when the width of the region, qL^, reaches 
Lx/2. This saturation occurs when the field reaches a 
critical field of 



Br 



ih/2e)iLx/dLi) = ih/2e){LxetVdp: 



(14) 



counterflow currents obey Eq. (11) and so jp{x) decays 



Beyond this critical field, all the phase domains in the 
sample take part in coherent tunneling with a constant 
magnitude for the tunneling current density jt, as de- 
picted in Fig. [2] (bottom). The counterflow current, 
therefore, rises linearly from zero from either edge, reach- 
ing a maximum at the center of the sample. This satu- 
rated current profile stays the same for fields higher than 
Be- Note that, whereas the critical field Bq in the clean 
limit is a microscopic parameter independent of sample 
geometry. Be is proportional to the length of the sample 
in a disordered system. 

We have tested the critical-state model by comparing 
it with numerical minimization of the energy [Eq. ([I])] 
for a one-dimensional lattice model. To do this we add 
a dissipative dynamical term —77 to the right-hand side 
of Eq. ([6| and numerically find the resulting steady- 
state. We begin by finding the steady-state with the 
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boundary condition drj/dx = 0\x=o,l^, corresponding to 
an isolated bilayer in the absence of an in-plane field. 
We then slowly increase the boundary condition so that 
dri/dx\x=a,L^ = —q, corresponding to applying the in- 
plane field. The counterflow supercurrent is then ob- 
tained from the resulting phase profile by adding the dia- 
magnetic term, as described above. The resulting current 
profiles, averaged over disorder realizations, are shown as 
the circles in Fig. [2j and are seen to agree closely with 
the critical-state model. 

We will now discuss the magnetic moment generated 
by the diamagnetic response of the bilayer. Integrating 
the current profile gives the moment Mn : 



LxLyd 




B 

2Ba 



B <B, 

B>Br. 



(15) 



We see that the moment no longer rises linearly with B 
with the full diamagnetic susceptibility xo of clean 
system. This is because, as seen in Fig. [2] the fraction of 
the sample with the full diamagnetic current jd is con- 
tinuously reduced as B (and hence q) increases. In fact, 
the total moment saturates at the critical field B^ to the 
value of 



Mil 



epsLx , 
-jfTjjLxLyd. 



(16) 



The analysis above applies if the field and supercurrents 
are sufficiently small that their presence does not mod- 
ify the pinning length L^, i.e., the supercurrents do not 
depin the vortices. This is the case if g < 1/^, the in- 
verse disorder correlation length, so that the additional 
phase winding at wavevector q in Eq. ([s]) does not affect 
the tunneling energy. For larger fields however, the tun- 
nel currents will oscillate within each correlation area, 
rapidly suppressing the net tunneling at larger scales. 
The current jp, which is jd at the boundary, will then be 
approximately uniform, and jcf and jt will be small at 
long wavelengths. Thus we expect the susceptibility to 
fall quickly for fields beyond 



B. 



c2 



h 

Ids, 



0.5T. 



(17) 



This can be regarded as a depinning field where the phase 
is no longer pinned by the random field but is deter- 
mined by the diamagnetic wavevector q. A correspond- 
ing suppressiorPsl of the current in interlayer tunneling 
experiments is expected, and observed, at such fields. 
Since Bc/Bc2 — Lx£,/2L'^ < 1 for experimental samples 
with Lx ^ 1 mm, the saturation at Be occurs before this 
suppression at Bc2- Note that in the lattice model dis- 
cussed above, the disorder is uncorrelated between sites, 
so that the numerical simulations do not capture the de- 
pinning effect at Bc2- 



IV. FLUX CREEP AND DISSIPATION AT LOW 
TEMPERATURES 

Let us now consider thermal fluctuations over the en- 
ergy barriers which pin the winding of the superfluid 
phase. This allows the supercurrents to relax to the 
equilibrium state with jcf = and jt — 0. The diamag- 
netic moments induced by an in-plane field will decay as 
phase slips enter the system. This gives rise to resistiv- 
ity in transport measurement^^. We will estimate the 
size of these effects, applying the conventional fiux-creep 
model of superconductorgSS. We neglect the possibility 
of vortex-glass stateP^^, which may further suppress the 
dissipation at low temperatures. 

For sufficiently low temperatures T, the dynamics of 
the merons controlling the diamagnetic moment will in- 
volve thermal fluctuations of the phase domains of size 
id discussed above. As in the Anderson-Kim theorjPH^ 
uncorrelated vortex motion should be irrelevant at low 
temperatures and small bias, because the vortex separa- 
tion ^ ^ Ld, so the energy barriers are larger. A thermal 
fluctuation in which the phase of a domain changes by 
approximately 27r corresponds to a phase slip moving a 
distance L^. This will occur at a frequency u!oe~^^^^^'''^ , 
where ujq is an attempt frequency and U{j) is an energy 
barrier. Such phase slip dynamics implies an interlayer 
voltage according to the Josephson relation V = W/e. 
For a typical domain size of Ld, this gives a layer- 
antisymmetric in-plane electric fleld 



bottom 



eid 



(18) 



At zero current, the typical energy barrier U{j) will be 
the energy of a domain, which is the stiffness ps irre- 
spective of siz^^H jji dimensions. The barrier will 
vanish at the current scale associated with the domain, 
jc — Id/Ld ~ eps/{hLd). The linear interpolation, 
U{j) = ps{l — j/jc), corresponds to the Anderson-Kim 
model. Inserting this form in Eq. (18), we obtain for the 



ohmic regime a sheet resistance for counterflow currents 
of 



Rs ~ {h^ujQ/e^kT)e-P^''''^ . [kT < p,) 



(19) 



The in-plane magnetic moment, and the long- wavelength 
counterflow supercurrent, relaxes in time t a^^ 



M|| (t)/M|i T=o « 1 - [kT/p,) Hioot), 



(20) 



where Mm 7-^0 is the T = moment from Eq. ( 15 ) 



It is possible to extend this argument into the dissipa- 
tive regime at high bias or temperature, where the barri- 
ers become irrelevant. The simplest assumption would be 
that the domain rotates every attempt time I/wq when 
the current density is j^.. This gives a resistivity in the 
flux-flow regime of 



9 ' 



(21) 
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so that we can rewrite the flux-creep form as 

Rs « Rle-P"''''^ . (22) 

However, details of the dissipation are likely to change be- 
tween these two regimes and other parameters may enter 
i?* , so that the prefactor in Eq. ( 22 ) will not be exactly 
the resistivity in the dissipative regime. Provided that 
R* is not fitted too far from the flux-creep regime, this is 
unimportant in practice because the exponential factor 
in general dominates in Eq. (22). Over a wider range 



of temperatures and currents, effects such as the activa- 
tion of quasiparticles may appear, leading to additional 
exponential factors in the dissipation. 

Reliable calculations of wq are very difficult, and ir- 
relevant in practice in the flux-creep form, Eq. (22 1. 



Nonetheless, an estimate may be obta ined b y considering 
the dynamical equation for the phasCj'i^Illl 



A 



2 sin(?7+6'o) 



e^Ps 



cq 



Pz 



v- 



Px 



(23) 



where pxx {pz) is the resistivity of in-plane (tunneling) 
quasiparticle currents, and c is the interlayer capacitance 
per unit area. This gives several candidates for the fre- 
quency scale at length scale L^. Since the resistivities are 
activated, we expect the first (incrtial) term in Eq. ( 23 ) 
to control the dynamics at sufhciently low temperatures: 
iuo = {e/hU){pJc)^/\ 



V. DISCUSSION 

We have examined the diamagnetic response of an iso- 
lated bilayer due to counterflow superfluidity. We ob- 
tained an in-plane diamagnetic moment which saturates 
at a critical value B^- Our theory also predicts a field 
scale Bc2 for the suppression of the in-plane diamagnetic 
response. There would be no in-plane susceptibility in a 
field-cooled bilayer, even below Bc2, allowing a separa- 
tion of the superfluid signal. We will now discuss briefly 
the magnitude and timescales of these effects using real- 
istic parameters for current experiments. 



From Eqs. (201 and (22 1, we see that the diamag- 



netic moment is long-lived, and the dissipation small, 
when kT Ps- In previous worlP^, we estimated that 
Ps « 20 mK in current experiments. This small value 
arises from the finite interlayer separation, and the re- 
duced area of the sample containing superfluid in the 
coherence network picture. Using this estimate as well 
as a similar area reduction for the capacitance, we find 
for the attempt frequency loq 300MHz. Since the su- 
perfiuid stiffness ps is around the lowest temperatures 
achieved, we expect from Eq. (20) that the moments re- 



lax rapidly. Similarly, counterfiow dissipation due to flux 
motion would be significant with i?* « 1 - 10 Ohm. 



Activated forms similar to Eq. ( 22 ) have previously 



been obtained for the residual counterfiow resistivity 
based on hopping of in-plane vortices (merons) within the 



coherence networlP^. In contrast, here we have obtained 
this form from the motion of phase-polarized domains, 
which corresponds to a large-scale collective motion of 
the merons. However, we expect the energy barriers for 
vortex motion in two dimensions to be the superfiuid stiff- 
ness up to logarithmic factors, so long as the possibility 
of a vortex glass is excluded. Thus this form is not de- 
pendent on the precise model of vortex motion. We note 
that experiments on hole bilayerdl^ see this activated be- 
havior in the counterflow resistance with a value of 20 
Ohms at 30 mK with an activated form that agrees with 
Eq. ( 19 1 if we use ps ~ 100 mK. Electron bilayer^ fit 



an approximately activated form with similar resistance 
values. A complete theory of dissipation and transport 
is beyond the scope of this work. Other mechanisms of 
dissipation will also play a role, such as the thermal ex- 
citation of vortices^ as well as fermionic quasiparticles 
with an energy gap of the order of ps ■ It suffices to point 
out that the contribution from fiux motion is not negli- 
gible in current samples. However the activated form of 
Eq. ( p2| suggests that a modest increase in p^, achiev- 
able by reducing the interlayer separation, would allow 
the study of the nearly dissipationless behavior of the 
fiux-creep regime. 

We now turn to the feasibility of measuring the dia- 
magnetic response of a bilayer at low temperatures. 
To do this, we compare the maximum moment My ^^ax 
[Eq. (16l] with the scale of the perpendicular moment, 
M± = LxLyeuJc/2, for an electron gas in the integer 
quantum Hall regime with cyclotron energy hujc ~ 50K 
in this system. Using our previous estimates Ps ~ 
20 mK, and Ld ~ 100 ^m, we obtain M|| „iax/-^-L = 
{ps/2njjJc){Lxd/L\) - L^/(2000m) when d = 28 nm. 
This smallness of the effect, in comparison to the con- 
ventional magnetic moment, makes this very challenging 
to measure using torque magnetometry. This can also be 
inferred by noting that in our theory the critical currents 
in tunneling experiments are the maximum diamagnetic 
currents, and the former are experimentally in the nano- 
ampere range. 

However, we can exploit our understanding of the dis- 
ordered isolated bilayer to consider how this diamagnetic 
response can be increased. In particular, we note that 
-^11, max per unit volume increases with and with the 
number of polarized domains in the sample. Note that 
the interlayer tunneling can be increased by many or- 
ders of magnitude compared with the current samples 
by reducing the tunnel barrier. Increasing the tunneling 
strength reduces the domain size L^. Also, the narrower 
barriers possible with stronger tunneling allow a larger 
Ps ■ Our approach holds up to the point where the domain 
size Ld reaches the disorder lengthscale, ^ ~ 100 nm, and 
the maximum field Be just reaches the depinning field 
Bc2- At this point, we obtain a significant magnetic mo- 
ment ,jjax/-^-L ^ L^l (1 mm). Thus the diamagnetism 
of the exciton superfluid should be evident in samples 
with strong interlayer tunneling. 

Although in the weak-tunneling samples the moments 
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are small, the pinning and dynamics of the in-plane flux 
can still be probed in transport experiments. One conse- 
quence of the pinning picture is that the critical current 
in a tunneling experiment should not be affected by in- 
planc fields smaller than Bc2, in contrast to the clean 
limit where the critical current is suppressed at fieldJ^ 
-Bq < Bc2- This is because in the pinning picture it is 
the gradients of the in-plane flux density which drive the 
flux through the disorder, and while the injected currents 
in the tunneling geometry impose such gradients, a uni- 
form in-plane field does not. It would also be possible to 
study the dynamics of the in-plane flux in tunneling ex- 
periments. An interlayer voltage at one end of the sample 
introduces in-plane flux at a given rate, and measuring 
the interlayer voltages and currents reveals the subse- 
quent motiorP^l of this flux. The logarithmic relaxation 
of the flux-creep regime can lead to hysteresis, which is 
seen in tunneling experimentJi^. 

An interesting extension of our work would be to con- 
sider samples in the Corbino disk geometry. In such a 
geometry, a radial magnetic field could induce circulat- 
ing counterflow supercurrents, similar to the induction 
of linear counterflow supercurrents by an in-plane field 
in the open geometry of Fig. [l] In the open geometry, 
tunneling is necessary to close the current loop, which we 
have seen gives rise to a maximum magnetic moment at 
the field Be- In a Corbino disk, tunneling is not needed to 
close a circulating current loop, so there may be no sat- 
uration effect at Be and the full magnetic response may 
persist up to Bc2- However, there will still be radial vari- 
ations in the densities of applied flux and diamagnetic 
current, and maintaining the diamagnetic current loop 
requires that these variations are pinned. Thus it seems 
likely that, even in this geometry, there is a maximum 
magnetic moment determined by the sample dimensions 
and a pinning length. The relevant pinning length, how- 
ever, may be determined by parameters other than the 
tunneling strength. 



VI. CONCLUSIONS 

In summary, we have considered the diamagnetic re- 
sponse of a bilayer exciton superfluid to an in-plane mag- 
netic fleld, in the presence of in-plane vortices nucleated 
by charge disorder. We argue that at low temperatures, 
changes in the in-plane magnetic field induce circulating 
diamagnetic currents and hence persistent diamagnetic 
moments. The maximum moments which can be induced 
are determined by the pinning of the in-plane fiux by the 
disorder, which involves a characteristic lengthscale re- 
lated to the tunneling. At finite temperatures, thermal 
motion of the in-plane fiux will lead to transport resis- 
tivities and cause the diamagnetic moments to decay. In 
current samples, we find that the maximum moments are 
small, but samples with stronger tunneling could allow 
persistent exciton supercurrents to be probed by torque 
magnetometry. Such experiments would be analogous to 
the torsional oscillator experiments that are the definitive 
measure of superfluid fraction in helium. 

Finally, we note that very strong tunneling, or very 
large stiffness, may even be able to prevent the dis- 
order nucleating vorticeJ^, recovering the clean limit 
in which the magnetic response is also expected to be 
measurable^*. In this case the phase diagram of the bi- 
layer would then closely parallel that of a type-II super- 
conductor, with both a mixed state (as described here) 
and a clean state, experimentally distinguishable by the 
presence of hysteresis. 
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